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Abstract 

The processes pp — > + X — » £ + £~j + X [l = e, fi) and pp — > + 
X — > vwy + X are calculated to 0(a s ) for general ZZ*y and Z77 couplings. 
The impact of 0(a s ) QCD corrections on the observability of ZZ'y and Z77 
couplings in Z7 production at the Tevatron and the Large Hadron Collider 
(LHC) is discussed. 
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I. INTRODUCTION 



Experiments at high energy hadron colliders provide an excellent opportunity to probe 
the interactions of the electroweak gauge bosons. Recent observation of electroweak gauge 
boson pair production at the Fermilab Tevatron have provided further confirmation of the 
electroweak Standard Model (SM) and have significantly tightened the constraints on pos- 
sible non-standard model self-interactions of the electroweak gauge bosons . For 
instance, production ]T],[2j can be used to study the WWj vertex, WZ production can 
be used to probe the WWZ vertex, and W + W~ production is sensitive to both the 
WW') and WWZ vertex functions. Furthermore, efforts have also been made to search 
for evidence of non-zero ZZj and Z77 couplings in Z(— > £ + £~)j (£ = e, /x) |5],[7||| and 
Z{— > vv)') production ]5||J. These couplings vanish in the SM at tree level. However, if 
new interactions beyond the SM are responsible for non-zero ZZ7 or Z77 couplings, then 
Z7 production could provide a clean signal for new physics. 

Previous studies of ZZ7 and Z77 couplings in hadronic Z7 production were based on 
leading order calculations [[H]]. Next-to-leading order calculations of pp — > Zf production 
in the SM have shown that the NLO corrections are largest at high values of the photon 
transverse momentum and high values of the Z7 invariant mass |IT] , [T2|| . These are the 
same regions of phase space where the effects of non-standard ZZ'y and Z77 couplings are 
most pronounced |10[ . It is therefore important to include NLO corrections when probing 
for evidence of non-standard ZZ'y and Z77 couplings in hadronic Z7 production. This 
paper presents a calculation of hadronic Z(— > £ + £~)^ and Z(— > Pu)^ production to (D(a s ), 
including the most general non-standard ZZj and Z77 couplings, and the decay of the Z 
boson in the narrow width approximation. Because of the larger Z — » vv branching ratio, 
the Z{— > vv)') cross section is about a factor 3 larger than the combined Z[— > e + e~)7 and 
Z{— > /i + /i _ )7 rates. This results in substantially better limits on anomalous ZZ'y and ^77 
couplings in the Z — > vv channel 

To perform our calculation, we use the Monte Carlo method for NLO calculations de- 
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scribed in Ref. In the calculation the SM is assumed to be valid apart from anomalies 



in the ZZ7 and Z77 vertices. In particular, we assume the couplings of W and Z bosons to 
quarks and leptons to be given by the SM, and that there are no non-standard couplings of 
the Z7 pair to two gluons fH| . The gluon fusion process, gg — > Z7, is small in the SM fl5 



and is not considered in this paper. Section II briefly reviews the method used to carry out 
the calculation, and describes how anomalous ZZ7 and Z77 couplings are incorporated. 

In Sec. Ill, we discuss how NLO QCD corrections influence the photon transverse momen- 
tum distribution, and derive sensitivity limits for non-standard ZZ'y and Z77 couplings at 
NLO for the Tevatron and LHC center-of-mass energies for various integrated luminosities. 
The photon transverse momentum distribution is the observable most sensitive to anoma- 
lous couplings ||10|| , and is used by CDF and D0 to extract information on the ZZ^f and 
Z77 vertices PJ^JUD- At the Tevatron, QCD corrections are modest, and slightly improve 
the sensitivity bounds. In contrast, at LHC energies, the inclusive NLO QCD corrections 
are quite large at high photon transverse momenta in the SM and reduce the sensitivity to 
anomalous ZZ'y and Z77 couplings somewhat. The large QCD corrections are caused by 
a log squared enhancement factor in the qg — > Z'-yq partonic cross section at high photon 
transverse momentum (pr), and the large quark-gluon luminosity at LHC energies. As in 
W-f @, WZ [0,11], and W+W~ production ppl, the effect of the QCD corrections at 



high pt at the LHC can be reduced by imposing a 0-jet requirement when searching for 
anomalous couplings. Finally, summary remarks are given in Sec. IV. 

II. OVERVIEW OF THE CALCULATION 



The calculation presented here generalizes the results of Refs. [JTTJ] and [12| to include gen 



eral (non-standard model) ZZj and Z77 couplings. The calculation employs a combination 
of analytic and Monte Carlo integration techniques; details of the method can be found in 
Ref. ||13|| . The leptonic Z boson decays are incorporated in the narrow width approximation. 
In this approximation, radiative Z decay diagrams, and graphs in which a virtual photon 
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decays into a charged lepton pair can be ignored. Radiative Z decays are of little interest 
when probing for non-standard couplings and can be suppressed by a suitable choice of cuts 
(see Sec. IIIB). Furthermore, in the narrow width approximation it is particularly easy to 



extend the NLO calculation of hadronic Z7 production in Ref. |TT[ to include the leptonic 



decay of the Z boson. The charged leptons are assumed to be massless in our calculation. 

A. Summary of 0(a s ) Z7 production including leptonic Z decay 

At lowest order in the SM, hadronic Z7 production proceeds via the Feynman diagrams 
shown in Fig. [l|. Non-standard ZZj and Z77 couplings contribute via the graphs shown in 
Fig. H]. At the leading-logarithm level, there are additional contributions to Z7 production 
which come from photon bremsstrahlung processes such as qg — » Zq followed by photon 
bremsstrahlung from the final state quark. Although the process qg —>■ Zq is formally 
of 0(aa s ), the photon fragmentation functions are of order a/a s JH); thus the photon 
bremsstrahlung process is of the same order as the Born process. 

The NLO calculation of Z7 production includes contributions from the square of the 
Born graphs, the interference between the Born graphs and the virtual one-loop diagrams, 
and the square of the real emission graphs. The basic idea of the method employed here 
is to isolate the soft and collinear singularities associated with the real emission subpro- 
cesses by partitioning phase space into soft, collinear, and finite regions. This is done by 
introducing theoretical soft and collinear cutoff parameters, 5 S and 5 C . Using dimensional 



regularization [22], the soft and collinear singularities are exposed as poles in e (the number 
of space-time dimensions is N = 4 — 2e with e a small number). The infrared singulari- 
ties from the soft and virtual contributions are then explicitly canceled while the collinear 
singularities are factorized and absorbed into the definition of the parton distribution func- 
tions. The remaining contributions are finite and can be evaluated in four dimensions. The 
Monte Carlo program thus generates n-body (for the Born and virtual contributions) and 
(n + l)-body (for the real emission contributions) final state events. The n- and (n + 1)- 
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body contributions both depend on the cutoff parameters S s and 5 C , however, when these 
contributions are added together to form a suitably inclusive observable, all dependence on 
the cutoff parameters cancels. The numerical results presented in this paper are insensitive 
to variations of the cutoff parameters. 

Except for the virtual contribution, the 0(a s ) corrections are all proportional to the 
Born cross section. It is easy to incorporate the leptonic Z decays into those terms which 
are proportional to the Born cross section; one simply replaces <i<5" Born (gg — > Zj) with 
da Born (qq — > Zj — > or da Born (qq — > Zj — ■> z/z/7) in the relevant formulae. When 

working at the amplitude level, the Z boson decay is trivial to implement; the Z boson 
polarization vector, e M (A;), is simply replaced by the corresponding Z — > £ + £~ or Z — > vv 
decay current in the amplitude. Details of the amplitude level calculations for the Born and 
real emission subpro cesses can be found in Ref. [^3[ . 



The only term in which it is more difficult to incorporate the Z boson decay is the 
virtual contribution. Rather than undertake the non-trivial task of recalculating the virtual 
correction term for the case of a leptonically decaying Z boson, we have instead opted to 
use the virtual correction for a real on-shell Z boson which we subsequently decay ignoring 
spin correlations. When spin correlations are ignored, the spin summed squared matrix 
element factorizes into separate production and decay squared matrix elements. Neglecting 
spin correlations slightly modifies the shapes of the angular distributions of the final state 
leptons, but does not alter the total cross section as long as no angular cuts (e.g., rapidity 
cuts) are imposed on the final state leptons. For realistic rapidity cuts, cross sections are 
changed by typically 10% when spin correlations are neglected. Since the size of the finite 
virtual correction is only about 2 — 4% the size of the Born cross section, the overall effect 
of neglecting the spin correlations in the finite virtual correction is expected to be negligible 
compared to the combined 10 — 20% uncertainty from the parton distribution functions, the 
choice of the factorization scale Q 2 , and higher order QCD corrections. 
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B. General ZZ^ and Z77 Couplings 



In qq — > Z7, the timelike virtual photon and/or Z boson couples to essentially massless 
fermions, which ensures that effectively <9 M V M = 0, V = 7, Z. This fact, together with gauge 
invariance of the on-shell photon, restricts the tensor structure of the Z'jV vertex sufficiently 
to allow only four free parameters. The most general non-standard ZZ'-f vertex function 
(see Fig. for notation) is P| 



m 



z 



rri z v ' 



+ hfe^q 2p (1) 

+ pa^ap , \ 

m| J 

where mz is the Z boson mass. The most general Z77 vertex function can be obtained from 
Eq. (^|) with the following replacements: 
p2 _ 2 p2 

^ — > and hf — »■ h] , % = 1, . . . , 4 . (2) 

mr z mr z 

Terms proportional to P M and q" have been omitted in Eq. (|IJ) since they do not contribute 
to the cross section. Without loss of generality, the overall ZZ'-f and Z77 coupling has been 
chosen to be 

9zz~f = 9 Zll = e , (3) 

where e is the charge of the proton. The overall factor (P 2 — qf) in Eq. ([I]) is a result of Bose 
symmetry, whereas the factor P 2 in the Z77 vertex function originates from electromagnetic 
gauge invariance. As a result, the Z77 vertex function vanishes identically if both photons 
are on-shell due to Yang's theorem p5| . 

All the anomalous couplings h\ (i = 1, . . . , 4, V = 7, Z) are C odd; h\ and h\ violate 
CP. hX and hY receive contributions only from operators of dimension > 8. Within the 
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standard model, at tree level, all the couplings h\ vanish. At the one loop level, only the 
CP conserving couplings h\ and h\ are nonzero |26[ . 



For simplicity, the Z boson mass mz has been chosen in Eq. ([I]) as the energy scale in the 
denominator of the overall factor and the terms proportional to /i^ 4 . For a different mass 
scale, M, all subsequent results can be obtained by scaling h\ z (^^ 4 ) by a factor M 2 /m| 
(M 4 /m|). 

Tree level unitarity restricts the ZZ^ and Z77 couplings uniquely to their standard 
model values at asymptotically high energies ||27| . This implies that the Z^V couplings 
hj have to be described by form factors hY(q 2 , q 2 , P 2 ) which vanish when q 2 , q 2 , or P 2 
becomes large. In Z7 production q\ = and q\ ~ m 2 z even when finite Z width effects 
are taken into account. However, large values of P 2 = s will be probed in future hadron 
collider experiments and the s dependence has to be included in order to avoid unphysical 
results that would violate unitarity. A detailed discussion of unitarity and form factors in 
nonstandard Z7 production can be found in Ref. flit)] . We will use generalized dipole form 
factors of the form 

h v 

hV(m 2 Zl 0, s) = ;° A2 . , (4) 

1 v z ' 1 + s A 2 FF ) n 



as advocated in Ref. [|10j. The subscript denotes the low energy value of the form factor. 
The mass scale App is the scale at which novel interactions, like multiple weak boson or 
resonance production, may appear. Unless stated otherwise, we use n = 3 (n = 4) for h\ z 
(h\ A ) and A FF = 500 GeV [A FF = 3 TeV] at the Tevatron [LHC]. 

At present, the most stringent direct bounds on anomalous Z^V couplings come from Z7 
production at the Tevatron collider, Z —>■ z/z/7 decays at LEP1, and 7+ invisible particles 
production at LEP2. From a search performed in the channels pp — > Z{— > £ + £~)7 and 
pp — > Z(—>- jyiy)^, the D0 Collaboration obtains from Run la data 

I /if, I < 0.78 (for hf = 0), \hf \ < 0.19 (for /if = 0), (5) 

and 



\tq \ < 0.81 (for hj = 0), \hj \ < 0.20 (for hj = 0), (6) 

at the 95% confidence level (CL). The limits obtained for h{ and virtually coincide with 
those found for h^ and h\ Q) respectively. The L3 Collaboration obtains a slightly better 
bound on h§ |[28|| . From e + e~ — > 7+ invisible particles at LEP2, the DELPHI collaboration 
finds a preliminary limit of \hj \ < 0.5 ||29|| . To derive these limits, the experiments assumed 
a form factor scale of App = 500 GeV with n = 3 (n = 4) for h\ (h\). 

It is straightforward to include the non-standard model couplings in the amplitude 
level calculations. We computed the qq — > virtual correction with the vertex func- 
tion of Eq. (|l|) in the 't Hooft-Veltman scheme p2],|3~0l using the computer algebra program 
FORM ||31|| . The resulting expression, however, is too lengthy to present here. 

Note that the non-standard ZZj and Z77 couplings do not destroy the renormalizability 
of QCD. Thus the infrared singularities from the soft and virtual contributions are explicitly 
canceled, and the collinear singularities are factorized and absorbed into the definition of 
the parton distribution functions, exactly as in the standard model case. 

The squared matrix element is bi-linear in the anomalous couplings. Due to the anti- 
symmetry of e fJ-a/3p , all terms proportional to h%, h\ , hf , and h\ hf (i = 1,2) vanish 
in the LO as well as in the NLO squared matrix elements, after the sum over the fermion 
helicities and the photon polarizations is performed. Terms proportional to h\ and h\ are 
present in the squared matrix elements. These terms are proportional to cos#*, where 9* is 
the scattering angle of the photon in the parton center of mass frame. They vanish after 
the integration over phase space is performed, unless rapidity cuts on the photon which are 
not symmetric with respect to rj = (where rj is the pseudorapidity) are chosen. 



III. PHENOMENOLOGICAL RESULTS 

We now discuss the phenomenological implications of 0(a s ) QCD corrections and general 
ZZ'y and Z77 couplings in Z7 production at the Tevatron (pp collisions at y/s = 1.8 TeV) 
and the LHC (pp collisions at yfs = 14 TeV). First, the input parameters, cuts, and the finite 
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energy resolution smearing used to simulate detector response are briefly described. We then 
discuss in detail the impact of 0(a s ) QCD corrections on the observability of non-standard 
ZZ'y and Z77 couplings in Z(— > £ + £~)7 and Z[— > vv)^ production at the Tevatron and 
LHC. As mentioned in the Introduction, we make no attempt to include the contributions 
from gluon fusion, gg — > Zj, which are formally of 0{a 2 ) 1 into our calculation. Gluon fusion 
contributes less than 0.2% (6%) to the total Z7 cross section at the Tevatron (LHC) |lq| . 



A. Input Parameters 

The numerical results presented here were obtained using the two-loop expression for 
a s . The QCD scale Aqcd is specified for four flavors of quarks by the choice of the parton 
distribution functions and is adjusted whenever a heavy quark threshold is crossed so that 
a s is a continuous function of Q 2 . The heavy quark masses were taken to be m& = 5 GeV 



and m t = 176 GeV 32.33 



The SM parameters used in the numerical simulations are mz = 91.187 GeV, mw = 
80.22 GeV, a(mw) = 1/128, and sin 2 = 1 — ( m w/ m z) 2 - These values are consistent 



with recent measurements at LEP, SLC, the CERN pp collider, and the Tevatron p4 
The soft and collinear cutoff parameters, discussed in Sec. IIA, are fixed to 5 S = 10~ 2 and 
5 C = 10~ 3 . The parton subprocesses have been summed over u, d, s, and c quarks. The Z 
boson leptonic branching ratio is taken to be B(Z — > e + e~) = 0.034 and the total width of 
the Z boson is Yz = 2.490 GeV. Except where otherwise stated, a sing le scale Q 2 = Mf r 
where M^ 7 is the invariant mass of the Zj pair, has been used for the renormalization scale 
fj 2 and the factorization scale M 2 . The NLO numerical results have been calculated in the 
modified Minimal Subtraction (MS) scheme . 



In order to get consistent NLO results it is necessary to use parton distribution functions 
which have been fit to next-to-leading order. The numerical simulations have been performed 
using the Martin-Roberts-Stirling (MRS) |5§ set A distributions (A 4 = 230 MeV) in the 
MS scheme. For convenience, the MRS set A distributions have also been used for the LO 
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calculations. 



B. Acceptance Cuts 

The cuts imposed in our numerical simulations are motivated by two factors: 1) the 
finite acceptance and resolution of the detector, and 2) the need to suppress radiative Z 
decays which result in the same final state as Z7 production but which are of little interest 
for the study of anomalous couplings in hadronic collisions. The finite acceptance of the 
detector is simulated by cuts on the four-vectors of the final state particles. This group of 
cuts includes requirements on the transverse momentum of the photon and charged leptons 
for Z(— > £ + £~)'~f, and on the missing transverse momentum, ]f> T , resulting from the non- 
observation of the neutrinos in Z{— ■> vv)^. Also included in this group are cuts on the 
pseudorapidity, i], of the photon and the charged leptons. In addition, the charged leptons 
and the photon are also required to be separated in the pseudorapidity-azimuthal-angle 
plane 



1 1/2 



A0 £7 ) +(A % ) . (7) 



Since we ignore photon radiation from the final state lepton line in our calculation, it is 
necessary to impose cuts which will efficiently suppress contributions from these diagrams. 
In radiative Z decays the lepton-photon separation sharply peaks at small values due to the 
collinear singularity associated with the diagrams in which the photon is radiated from the 
final state lepton line. In the following we shall therefore impose a large separation cut of 
AR(£, 7) > 0.7. Contributions from Z — > £ + £~ / y can be further reduced by an invariant 
mass cut on the ££^ system of M(££j) > 100 GeV. 

At leading order, Z'j events are produced not only by the Born subprocess qq — > Z^ 
but also by the photon bremsstrahlung process which proceeds via subprocesses such as 
qg — > Zq followed by photon bremsstrahlung from the final state quark. As demonstrated in 
Ref. ||37[| , the bremsstrahlung process is significant at LHC energies. However, this process 
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does not involve the Z^V vertices and is thus a background in a search for anomalous ZZ^i 
and Z77 couplings. Fortunately, the photon bremsstrahlung events can be suppressed by 
requiring the photon to be isolated j37|. A photon isolation cut typically requires the sum 
of the hadronic energy -Ehad in a cone of size Rq about the direction of the photon to be less 
than a fraction eh of the photon energy E 7 , i.e., 



A_R<_R 



(8) 



with AR = [(Acj)) 2 + (Ar/) 2 ] 1 / 2 . To suppress the photon bremsstrahlung background, a 
photon isolation cut with eh = 0.15 and Rq = 0.7 |JS] will be applied in the numerical results 
presented in this section. For this value of eh, the photon bremsstrahlung background is less 
than a few per cent of the Born Z7 signal rate. The complete set of cuts for Z{— > £ + £~) , ~f 
is summarized in the following table. 

Z 7 -> 

Tevatron LHC 



p T ( 7 ) > 10 GeV 
p T (£) > 15 GeV 
1^(7)1 < 3.0 
W)\ <3.5 
AR(£,-f) > 0.7 
M(Uj) > 100 GeV 

EAK<0.7 -^had < 0.15 .Ely 



p T (j) > 100 GeV 
p T {£) > 20 GeV 
1*7(7)1 < 3.0 

\n(£)\ < 3.0 

AR(£,j) > 0.7 
M(%) > 100 GeV 

Sa,R<0.7 -^had < 0.15 -Ely 



Since our calculation is carried out in the narrow width approximation for the Z boson, no 
explicit cut on the di-lepton invariant mass is imposed. 

If the Z boson decays into a pair of neutrinos, the experimental signal is pp — > 7^ T , with 
the missing transverse momentum, $ T , resulting from the nonobservation of the neutrino 
pair. For Z7 — > ^ T 7 at the Tevatron we use the same transverse momentum and pseudo- 
rapidity cuts as the D0 Collaboration in their Z[— > z/z/) 7 analysis |J. The following table 
summarizes the cuts imposed for both the Tevatron and LHC analysis. 
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Tevatron LHC 



p T (j) > 40 GeV p T (j) > 100 GeV 

i> T > 40 GeV i> T > 100 GeV 

no jet with pr(j) > 15 GeV and \r](j)\ < 2.5 no jet with pt{j) > 50 GeV and \rj(j)\ < 3 
17/(7)1 < 2.5 |77( 7 )| < 3.0 

The high pril) and ^ T cuts, when combined with the jet veto, strongly reduce the back- 
ground from W — > ev events where the electron is misidentified as a photon at the Tevatron. 
These cuts also eliminate backgrounds from jj production with the jet rapidity outside the 
range covered by the detector and thus "faking" missing transverse momentum, and jj 
production where one of the jets is misidentified as a photon, while the other disappears 
through the beam hole. The large p T {^) and f T cuts at LHC energies are chosen to reduce 
potentially dangerous backgrounds from Z + 1 jet production, where the jet is misidentified 
as a photon, and from processes where particles outside the rapidity range covered by the 



detector contribute to the missing transverse momentum. Present studies [39,40] indicate 



that these backgrounds are under control for p T (j) > 100 GeV and f T > 100 GeV. 

C. Finite Energy Resolution Effects 

Uncertainties in the energy measurements of the charged leptons and photon in the de- 
tector are simulated in the calculation by Gaussian smearing of the particle four-momentum 
vector with standard deviation a. For distributions which require a jet definition, e.g., the 
Z'-f + 1 jet exclusive cross section, the jet four-momentum vector is also smeared. The 
standard deviation a depends on the particle type and the detector. The numerical results 
presented here for the Tevatron and LHC center of mass energies were obtained using a 
values based on the CDF [JjJ and ATLAS [|lj specifications, respectively. For Z{— > Pi/) 7 



production at the Tevatron, the photon and f T vectors were smeared using the D0 resolu- 
tions given in Ref. |§. 
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D. 0(a s ) Corrections and Anomalous ZZ^ and Z77 Couplings in Z(— > £ + £ )j 

Production 

Non-standard ZZ7 and Z77 couplings have a significant effect on many distributions 



in Zj production. The photon transverse momentum distribution was found [10| to be the 
distribution most sensitive to anomalous couplings. We will therefore concentrate on it in 
the following presentation. 

The LO and NLO SM photon transverse momentum distributions for pp — > Z7 + X — ► 
£ + £~^ + X production at Tevatron and LHC center of mass energies are shown in Fig. 
Here, and in all subsequent results shown for Z — > £ + £~ decays, we sum over electron and 
muon final states. The NLO corrections grow with the photon transverse momentum and 
with the center of mass energy. At the LHC, for example, the QCD corrections increase 
the SM cross section by about a factor 2.2 at Pt(j) — 1 TeV, whereas the enhancement is 
only a factor 1.4 at Pt(j) = 100 GeV. The large QCD corrections at high values of Pt(i) 
are caused by a collinear enhancement factor, log 2 (pxij) / mz) , in the qg — > Z'jq partonic 
cross section for photon transverse momenta much larger than rriz, Pr{l) ^> mzi an d the 
large qg luminosity at LHC energies. The large corrections arise from the kinematical region 
where the photon is produced at large p T and recoils against the quark, which radiates a 
soft Z boson which is almost collinear to the quark, and thus is similar in nature to the 
enhancement of QCD corrections observed at large vector boson transverse momenta in 
Wj, WZ, and W + W~ production |l6| , [T8| , pO| . The effect, however, is less pronounced than 
in these processes. In Wj and WZ production, the SM Born cross section is suppressed 



due to the appearance of an exact or approximate radiation zero while there is 

no radiation zero in the Z7 case. In W + W~ production, the strong correlation of the W 
helicities in the SM, together with the effect of kinematic cuts, is responsible for the larger 
effect of the QCD corrections [19| . 



The effects of non-standard ZZ7 couplings on the photon transverse momentum distri- 
bution in pp — > Zj + X — » £ + £~7 + X at the Tevatron center of mass energy are shown in 
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Fig. [5| The LO and NLO results are shown in Fig. [|a and Fig. |5|b, respectively. Results are 
displayed for the SM and two sets of anomalous couplings, (h§ = 1.0, hf = hf = /if = 0, 
SM Z77 couplings) and (/if = 0.05, /if = hf = /if = 0, SM Z77 couplings). For simplic- 
ity, only one coupling at a time is allowed to differ from its SM value. In order to clearly 
display the effect of the anomalous couplings, we have chosen rather large values for /if and 
hf , here as well as for the LHC (see below). However, S"-matrix unitarity is respected for 
the chosen values of the anomalous couplings, the power of the form factor, and the form 
factor scale. The 0(a s ) corrections in the presence of anomalous couplings at the Tevatron 
energy are approximately 20 - 40%, as in the SM. 

For equal coupling strengths, the numerical results obtained for the Z77 couplings h1 Q 
and hj are about 20% below those obtained for hf and hf in the region where anomalous 
coupling effects dominate over the SM cross section. Results for the CP-violating couplings 
^12 = 7) are identical to those obtained for the same values of h% 4 . Since terms linear 
in the anomalous couplings vanish in the differential cross sections, results are insensitive to 
the sign of the anomalous couplings if only one coupling at a time is allowed to differ from 
its SM value. 

The Pt{i) distribution for Z(— > l + l~~) r y production at the LHC is shown in Fig. ^. At 
leading order, the sensitivity of the photon transverse momentum distribution to anomalous 
ZZ^ couplings is significantly more pronounced than at the Tevatron. In the presence of 
anomalous couplings, the higher order QCD corrections are considerably smaller than in the 
SM. For large values of Pt{i), when anomalous couplings dominate, the 0(a s ) corrections 
are typically between 20% and 40%. In the same region, QCD corrections enhance the SM 
cross section by about a factor 2.2. At next-to-leading order, the sensitivity of the photon 
transverse momentum spectrum to anomalous couplings thus is somewhat reduced at the 
LHC. The logarithmic factor causing the cross section enhancement at high values of Pt{i) 
in the SM originates from the collinear region. The Feynman diagrams contributing in 
this region do not involve the Z^V vertices. The logarithmic enhancement factor therefore 
does not affect the anomalous contributions to the matrix elements. Because /if receives 
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contributions only from operators with dimension > 8, terms in the helicity amplitudes 
proportional to it grow like (y/s/mz) 5 - Deviations originating from /if, therefore, start at 
higher invariant masses and rise much faster than contributions from couplings such as /if 
which correspond to dimension 6 operators. 

The effect of the QCD corrections is shown in more detail in Fig. [7], where we display 
the ratio of the NLO and LO differential cross sections for the transverse momentum of 
the photon. At the Tevatron, the NLO to LO cross section ratio slowly rises from 1.2 at 
Pt{i) = 10 GeV to about 1.5 at pr{l) = 400 GeV. Since we used a rather small form 
factor scale of App = 500 GeV for the Tevatron, the effect of the anomalous couplings is 
suppressed at high transverse momenta. As a result, the NLO to LO cross section ratio for 
non- vanishing anomalous couplings is very similar to that obtained in the SM. At the LHC, 
we use App = 3 TeV and the cross section ratio gradually decreases with pr{l) from pa 1.35 
to 1.2 if anomalous couplings are present. In contrast, the SM NLO to LO cross section 
ratio increases from pa 1.35 at Pr(l) = 100 GeV to 2.2 at Pt(i) = 1 TeV. 

From the picture outlined above, one expects that, at next-to-leading order, a large 
fraction of the Zj events with large photon transverse momentum will contain a high px 
jet at the LHC. At the Tevatron, on the other hand, the jet activity of Z7 events at high 
Pt{i) should be substantially reduced. This fact is illustrated in Fig. || which shows the 
decomposition of the inclusive SM NLO Pt(j) differential cross section into NLO 0-jet and 



LO 1-jet exclusive cross sections 0]. For comparison, the photon transverse momentum 
distribution in the Born approximation is also shown in the figure. Here, a jet is defined as 
a quark or gluon with 

PtU) > 10 GeV and \r]{j)\ < 2.5 (9) 

at the Tevatron, and 

PtU) > 50 GeV and \rj(j)\ < 3 (10) 

at the LHC. The sum of the NLO 0-jet and the LO 1-jet exclusive cross section is equal 
to the inclusive NLO cross section. The NLO exclusive Z7 + jet and the LO exclusive 
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Z^i + 1 jet cross sections depend explicitly on the jet definition, however, the inclusive NLO 
cross section is independent of the jet definition. 

Present LHC studies [[^,[4(].[44|] and projections to Tevatron energies suggest that jets 



fulfilling the criteria of Eqs. (|9|) and (|1(]) can be identified in Z7 + X events at the Teva- 
tron [35] and LHC [3JJ] for luminosities up to 10 33 cm -2 s _1 and 10 34 cm -2 s -1 , respectively. 
It should be noted, however, that for theoretical reasons, the jet transverse momentum 
threshold can not be made arbitrarily small in our calculation. For transverse momenta 
below 5 GeV (20 GeV) at the Tevatron (LHC), soft gluon resummation effects are expected 



to significantly change the shape of the jet p T distribution \^7j . For the jet definitions dis- 
cussed above, these effects are expected to be unimportant and therefore are ignored in our 
calculation. 

Figure |8| shows that, at the Tevatron, the 1-jet cross section is always considerably smaller 
than the NLO 0-jet rate. At the LHC, on the other hand, the 1-jet cross section is larger than 
the LO cross section for Pt(i) > 400 GeV. The effect of the QCD corrections can be reduced 
by vetoing hard jets in the central rapidity region, i.e., by imposing a "zero jet" requirement 
and considering Z7 + jet production only. The NLO 0-jet and Born differential cross 
sections deviate by 30% at most in the Pt region shown. The photon transverse momentum 
distribution for NLO Z7 + jet production is shown in Fig. || The 0-jet requirement is seen 
to restore the sensitivity to anomalous couplings lost in the inclusive NLO cross section at 
the LHC. It has little effect at the Tevatron. 

As mentioned in Sec. IIIA, all our results are obtained for a scale of Q 2 = M§ . The 
Born cross section for Z7 production depends significantly on the choice of Q, which enters 
through the scale-dependence of the parton distribution functions. At the NLO level, the 
Q-dependence enters not only via the parton distribution functions, but also through the 
running coupling a s (Q 2 ) and the explicit factorization scale-dependence in the order a s (Q 2 ) 
correction terms. Similar to the situation encountered in W7, WZ, and W + W~ production 
in hadronic collisions |TB|,[T^,[TP[ , we find that the NLO Z7 + jet exclusive cross section is 



almost independent of the scale Q. Here, the scale-dependence of the parton distribution 
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functions is compensated by that of a s (Q 2 ) and the explicit factorization scale dependence 
in the correction terms. The Q-dependence of the inclusive NLO cross section is larger than 
that of the NLO 0-jet cross section; it is dominated by the 1-jet exclusive component which 
is calculated only to lowest order and thus exhibits a considerable scale-dependence. 



E. 0{a s ) Corrections and Anomalous ZZ7 and Z77 Couplings in Z{— > 1/^)7 Production 

If the Z boson produced in qq — > Zj decays into neutrinos, the signal consists of a high 
Pt photon accompanied by a large amount of missing transverse momentum, j> T . Due to 
the larger Z — > vv branching ratio, the r )j> T + X differential cross section is about a factor 3 
larger than that for qq — > e + e~7 + X and qq — > /! + /i~7 + X combined. This results in limits 
on the anomalous ZZ'j and Z77 couplings which are presently about a factor 2 better than 
those obtained from the Z{— > £ + £~)-f analysis ||. 

The NLO photon transverse momentum distribution for Z{— > vv)^ production at the 



Tevatron is shown in Fig. |T0|a. Here we have imposed the cuts of the DO Z{— > vv)^ analysis 
(see Sec. IIIB and Ref. [||]). Since jets with a transverse energy larger than 15 GeV are 
excluded in the experimental analysis, only the 0-jet differential cross section is shown. 
The effect of non-standard Z^V couplings is very similar to that observed in Z[— > l + £~) r ) 
production. The impact of the QCD corrections on the differential cross section is shown in 
Fig. |T0|b, where we display the ratio of the NLO and LO differential cross sections for the 
transverse momentum of the photon. Due to the jet veto cut imposed, the QCD corrections 
are small over a wide range of photon transverse momenta. For Pt(i) < 200 GeV, the cross 
section ratio is almost constant. It rises slowly for larger values of Pr{l)- The NLO to LO 
cross section ratio in the SM and in the presence of anomalous couplings are very similar. 
Figure |T0|b demonstrates that, for data samples containing only a few Z(— > ui/)j events 



with pr(7) > 200 GeV and for the jet veto cut used, the LO calculation is an adequate 
approximation. However, increasing or decreasing the jet transverse momentum threshold 
will change the size of the QCD corrections, and thus increase the deviation of the NLO 
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calculation from the LO result. The calculation of Ref. , where the effect of the NLO QCD 



corrections is approximated by a constant fc-factor of k — 1 + 8a s /9ir « 1.34, overestimates 
the cross section by about 30% over a wide range of photon transverse momenta. This 
calculation has been used in the D0 analysis to compare data with the SM prediction, and 
to extract limits for the ZZ^ and Z77 couplings. Since it overestimates the cross section, 
the limits obtained are slightly better than those extracted using the full NLO calculation 
(see Sec. IIIF). 

In Fig. [Tl| a we show the NLO photon transverse momentum distribution for Z{— > vv)^ 
production at the LHC, imposing a pr(j) < 50 GeV jet veto cut. Figure |ll]b displays the 
ratio of the NLO and LO differential cross sections for the transverse momentum of the 
photon. For the pr(j) threshold chosen, NLO QCD corrections reduce the cross section by 
up to 20%. In contrast to the situation encountered at the Tevatron, the cross section ratio 
slowly falls with pril)- 

F. Sensitivity Limits 

We now study the impact that 0{ot s ) QCD corrections to Z7 production have on the 
sensitivity limits for hj at the Tevatron and LHC. For the Tevatron we consider integrated 
luminosities of 1 fb _1 , as envisioned for the Main Injector era, and 10 fb _1 (TeV33) which 
could be achieved through additional upgrades of the Tevatron accelerator complex |^5| . In 
the case of the LHC we use fCdt = 10 fb -1 and 100 fb _1 To extract limits in the 

Z{— > £ + £~) , ~f case, we sum over electron and muon final states. Interference effects between 
different Z^jV (V = Z, 7) couplings are fully incorporated in our analysis. 



assum- 



To derive 95% CL limits we use the Pt(i) distribution and perform a \ test H3 
ing that no deviations from the SM predictions are observed in the experiments considered. 
We include the cuts summarized in Sec. IIIB. For Z{— > £ + £~) / y production, we use the jet 
definitions of Eqs. (|9|) and (|T0"D. Unless explicitly stated otherwise, a form factor as given 
in Eq. (^) is used with n = 3 for hY 3 , and n = 4 for /i^ 4 . Furthermore, the form factor 
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scale Aff is taken to be 0.5 TeV (3.0 TeV) for Tevatron (LHC) simulations. The Pt{i) 
distribution is split into a certain number of bins. The number of bins and the bin width 
depend on the center of mass energy and the integrated luminosity. In each bin the Poisson 
statistics are approximated by a Gaussian distribution. In order to achieve a sizable count- 
ing rate in each bin, all events above a certain threshold are collected in a single bin. This 
procedure guarantees that a high statistical significance cannot arise from a single event 
at large transverse momentum, where the SM predicts much less than one event. In order 
to derive realistic limits we allow for a normalization uncertainty of 50% in the SM cross 
section. For the cuts we impose, background contributions other than SM Z'j production 



are small |H| and are ignored in our derivation of sensitivity limits. 

The calculation of sensitivity bounds is facilitated by the observation that the CP con- 
serving couplings /134 and the CP violating couplings h\ 2 do not interfere. Furthermore, 
cross sections and sensitivities are nearly identical for equal values of hi 02 o and h^ 040 . In 
the following we shall therefore concentrate on h^ 4 . In each bin, i, the cross section is a 
bilinear function of the anomalous couplings: 

o i = o i (SM) + £ (af hV + a?hJ)+ £ £ b% v ' h] h v k ' . (11) 

V=7, Z V,V'=-y, Z j,k=3A 

Here, a l (SM) is the SM cross section, and af and b^. v are constants. Since the interference 
terms between the SM and the anomalous contributions to the helicity amplitudes vanish 
for h\ 4 (see Sec. IIB), 

af = af = 0. (12) 

These constraints are taken into account in our calculation of sensitivity bounds. 

For hX 2 , an expression for the cross section similar to that of Eq. ( |TTD can be derived. In 



this case, however, the coefficients of the terms linear in h\ 2 only vanish if the phase space 
integration over the scattering angle of the photon, 9, is symmetric in cos# (see Sec. IIB). 

Our results are summarized in Figs. [TJ| - |T5| and Tables | - JTJ. Figure |T2| shows 95% CL 
contours in the /if — hj , /if — hJ Q , and the hf — hj plane for Z{— > production 
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and an integrated luminosity of 1 fb 1 at the Tevatron. Results for the Kl — hf plane are 
very similar to those displayed in Fig. [T2| b for h§ versus hj , and are therefore not shown. 



In each figure, only those couplings which are plotted against each other are assumed to 
be different from their zero SM values. As noted in Ref. flOf, ZZ'-f and Z77 couplings 



interfere little at LO. Figure ^demonstrates that the NLO QCD corrections do not change 
this behaviour. This statement also applies to different integrated luminosities and to Z{— > 
£ + £~) , y production at the LHC. Due to the larger cross section, the sensitivity bounds 
obtained from the inclusive NLO differential cross section are about 5% better than those 
derived using the LO calculation. The limits extracted from the NLO Z7 + O jet cross section 
are almost identical to those found using the LO calculation. 

Larger correlations are encountered between h\ and h\ (see Ref. |10|). The impact of 



the 0(a s ) QCD corrections on the limits in the hf — hf plane at the Tevatron is shown in 
Fig. [13] for integrated luminosities of 1 fb -1 and 10 fb -1 . The inclusive NLO QCD corrections 
improve the sensitivity bounds by up to 9% (6%) for 1 fb -1 (10 fb -1 ) whereas the limits 
obtained analyzing the Z7 + jet channel are very similar to those found using the LO 
calculation. Similar results are obtained for hi and hj. The sensitivity bounds for the Z77 
couplings are a few per cent weaker than those found for the corresponding ZZj couplings, 
and QCD corrections have a smaller effect than in the case of the ZZ7 couplings. Table | 
summarizes the 95% CL sensitivity limits, including all correlations, at leading order and 
next-to-leading order for hj and h\ Q {V = Z, 7) for the process pp — ► Z^ + X — > £ + £~j + X 
at the Tevatron with J Cdt = 1 fb -1 . 

Figure [14] displays how the 95% CL contour limits for pp — > Z7 + X — > ^> T j + X at 
the Tevatron and an integrated luminosity of 1 fb -1 are affected by NLO QCD corrections. 
Since a jet veto cut of Pt{j) < 15 GeV is imposed, we only show the LO (solid line) and NLO 
0-jet (dotted line) contour limits. As we have demonstrated in Sec. HIE, the NLO QCD 
corrections to Z(—> vv)^ production at the Tevatron are small over a large range of photon 
transverse momenta due to the 0-jet requirement. The sensitivities achievable for hf Q and 
hf in pp — > Z^i + X — > ^ T 7 + X at the Tevatron for / Cdt = 1 fb -1 and / Cdt = 10 fb -1 are 
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listed in Table ||. Besides the LO and NLO 0-jet bounds, we also show the limits obtained 
using the calculation of Ref. JED] where the NLO corrections are approximated by a constant 
fc-factor given by k = 1 + 8a s /97r « 1.34. This calculation has been used to extract bounds 
on Z'yV couplings from Z{— > vv)^ production at the Tevatron ||. Since a constant fc-factor 
overestimates the cross section for the jet veto cut imposed in the current experimental 
analysis of the Z(— > vv)^ channel, the sensitivity limits obtained are 10 — 13% better than 
those found using the full NLO 0-jet calculation, depending on the integrated luminosity. 
Results similar to those shown in Fig. |14] and Table [ij are obtained for /13 4 . 



The 95% CL limit contours in the h^ — hf plane for Z{— > £ + £ )j production at the LHC 
are shown in Fig. 15, assuming an integrated luminosity of 100 fb _1 . Table |T| summarizes 



the LO and NLO sensitivity bounds for pp — > Z7 + X — > £ + £~7 + X and pp — > + X — > 
^ T 7 + X at y/s = 14 TeV with JCdt = 10 fb _1 and 100 fb _1 . At LHC energies, the inclusive 
0(a s ) QCD corrections in the SM considerably change the shape of the Pt(j) distribution 
(see Fig. ||b). As a result, the inclusive NLO QCD corrections reduce the sensitivity to 
anomalous couplings by 7 - 10%. As the integrated luminosity increases, larger transverse 
momenta become accessible. The difference between the LO and NLO sensitivity bounds 
for 100 fb" 1 therefore is slightly larger than for 10 fb _1 . In Sec. HID we have demonstrated 
that the size of the 0(a s ) QCD corrections at the LHC in the high pril) region can be 
reduced by vetoing hard jets in the central rapidity region. The sensitivity bounds obtained 
for the Z7 + jet channel are about 5% better than those found for the inclusive NLO case. 
However, because the NLO 0-jet cross section is smaller than the LO cross section for the jet 
definition we use (see Fig. |8]b), the limits obtained in the NLO Z7 + jet case are slightly 
worse than those extracted from the LO cross section. 

As we have mentioned in Sec. HID, the NLO Z7 + jet differential cross section is more 
stable to variations of the factorization scale Q 2 than the LO and inclusive NLO Z'j + X 
cross sections. The systematic errors which originate from the choice of Q 2 will thus be 
smaller for bounds derived from the NLO Z7 + jet than those obtained from the inclusive 
NLO Z^i + X or the LO cross section. 
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As discussed in Ref. ||10|| , the limits which can be achieved are sensitive to the form and 
the scale of the form factor. For example, doubling the form factor scale to App = 1 TeV 
at the Tevatron improves the bounds by almost a factor 3. The dependence of the limits 
on App can be understood easily from Figs. [5] and || The improvement in sensitivity with 
increasing A FF is due to the additional events at large Pt(i) which are suppressed by the 
form factor if the scale App has a smaller value. To a lesser degree, the bounds also depend 
on the power n in the form factor. Reducing n allows for additional high Pt(t) events and 
therefore leads to a somewhat increased sensitivity to the low energy values of the anomalous 
couplings. It should be noted, however, that n must be larger than 3/2 (5/2) for h\ z (^4) 



in order to preserve S'-matrix unitarity ||T0 |. 



The bounds derived in this section are quite conservative. Using more powerful statistical 
tools than the simple x 2 test we performed can lead to considerably improved limits f|8| . 



The effect of the NLO QCD corrections on the sensitivity bounds, however, does not depend 
on the technique used to extract them. 

IV. SUMMARY 

Hadronic Z7 production can be used to probe for non-standard self interactions of 
the photon and Z boson. The experimental limits for non-standard ZZ'-f and Z77 cou- 
plings ||[7],[8]J| so far have been based on leading order calculations ffTOfl . In this paper we 



have presented an 0{a s ) calculation of the reactions pp — > Z7 + X — > £ + £~7 + X and 
PJ? — > Z^i + X — > ^ T 7 + X for general ZZ^i and Z77 couplings. The leptonic decay of the 
Z boson has been included in the narrow width approximation in our calculation. Decay 
spin correlations are correctly taken into account, except in the finite virtual contribution. 
The finite virtual correction term contributes only at the few per cent level to the total NLO 
cross section, thus decay spin correlations can be safely ignored here. Photon radiation from 
the final state lepton line is not taken into account; effects from Z —>■ £ + £~ / y decays can 
easily be suppressed by imposing a ££7 invariant mass cut and a cut on the lepton photon 
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separation. 

The Pt(t) differential cross section is very sensitive to nonstandard Z7V (V = Z, 7) 
couplings. QCD corrections change the shape of this distribution. The shape change is 
due to a logarithmic enhancement factor in the qg and qg real emission subprocesses which 
appears in the high Pt{i) region of phase space where the photon is balanced by a high 
Pt quark which radiates a soft Z boson. The logarithmic enhancement factor combined 
with the large gluon density at high center of mass energies make the 0(a s ) corrections 
quite large for Pt{i) 3> rriz- Since the Feynman diagrams responsible for the enhancement 
at large pri'j) do not involve any Z'-fV couplings, inclusive 0(a s ) QCD corrections to 
production tend to reduce the sensitivity to anomalous couplings. 

At the Tevatron, Z7 production proceeds mainly via qq annihilation. Here the main 
effect of the QCD corrections is an increase of the cross section by about 20 - 25%. The 
sensitivity limits derived from the inclusive NLO Z7 + X cross section at the Tevatron 
are thus up to 10% better than those obtained from the LO cross section. If a jet veto 
is imposed, the limits are almost identical to those found using the LO calculation. In its 
pp — > Z^f — > ^ T 7 analysis, the D0 Collaboration imposes &Pt(j) < 15 GeV requirement, and 
uses the calculation of Ref. JIT]] where the effect of the NLO QCD corrections is approximated 
by a simple fc-factor to extract sensitivity limits. We found that the bounds obtained from 
the full NLO Z{— > vv)j + jet calculation are 10 — 13% weaker than those derived using 
the calculation used in the experimental analysis. 

At the LHC, qg fusion significantly contributes to Z7 production and the change in the 
slope of the Prin) distribution caused by the NLO QCD corrections is quite pronounced. As 
a result, the limits on ZZ7 and Z77 couplings extracted from the inclusive NLO Z7 + X 
cross section are up to 10% weaker than those extracted using the LO calculation. The size 
of the QCD corrections at large photon transverse momenta can be reduced considerably, 
and a fraction of the sensitivity to Z7 couplings which was lost at the LHC may be regained 
by imposing a jet veto. The improvement, however, is moderate. 

Although a jet veto does not have a large effect on the sensitivity bounds at the Tevatron 
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or the LHC, extracting limits from the Z7 + jet channel has the advantage of a reduced 
uncertainty from the variation of the factorization scale Q 2 ; the dependence of the NLO 
Z'j + jet cross section on Q 2 is significantly smaller than that of the inclusive NLO and 
the LO Z'-f cross section. 

The effect of QCD corrections on the sensitivity limits for anomalous gauge boson cou- 
plings in Z7 production is significantly smaller than for Wj, WZ, and W + W~ production. 
In and WZ production, the SM Born cross section is suppressed due to the appearance 
of an exact or approximate radiation zero |)42ft3| . In W + W~ production, the strong correla- 
tion of the W helicities in the SM, together with the effect of kinematic cuts, is responsible 



for the larger effect of the QCD corrections [19 
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TABLES 

TABLE I. Sensitivities achievable at the 95% confidence level for the anomalous Z^V couplings 
h% and h\ (V = Z, 7) in pp -> Z7 + X + X, t = e, /x, at the Tevatron (y/s = 1.8 TeV) 

with J Cdt = 1 fb _1 . The limits for each coupling apply for arbitrary values of the other couplings 
listed in this table. The CP violating couplings hY 2 are assumed to take their SM values. For the 
form factor we use the form of Eq. (||) with n = 3 (n = 4) for (hY) and Aff = 0.5 TeV. The 
cuts summarized in Sec. IIIB are imposed. For the jet definition, we have used Eq. (^). 



coupling 


LO 


NLO incl. 


NLO 0-jet 


\h z 1 
l n 30l 

\h z 1 


0.62 
0.136 


0.58 
0.124 


0.61 
0.130 


l n 30l 

H \ 


0.65 
0.141 


0.64 
0.138 


0.68 
0.148 
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TABLE II. Sensitivities achievable at the 95% confidence level for the anomalous ZZ7 couplings 
h§ and hf in pp — ► Z7 + X — > ^ T 7 + ^ at the Tevatron (y^i = 1.8 TeV) for a) / Cdt = 1 fb _1 , 
and b) / Cdt = 10 fb _1 . Shown are the limits obtained from the LO calculation, the full NLO 0-jet 



differential cross section, and the calculation of Ref. [IV]. The limits for each coupling apply for 
arbitrary values of the other coupling listed in this table. The CP violating couplings hf 2 and all 
Z77 couplings are assumed to vanish. For the form factor we use the form of Eq. (Q) with n = 3 
(n = 4) for h§ (hf) and Aff = 0.5 TeV. The cuts summarized in Sec. IIIB are imposed. For the 
jet definition, we have used Eq. (0). 







a) JCdt = 1 fb" 1 




coupling 


LO 


NLO 0-jet 


NLO appr. 


\h z 1 
r<-3ol 


0.55 


0.53 


0.46 


\h z 1 


0.108 


0.104 


0.091 






b) JCdt = 10 ft" 1 




coupling 


LO 


NLO 0-jet 


NLO appr. 


1*1,1 


0.30 


0.29 


0.26 


\h z 1 

l n 40l 


0.055 


0.053 


0.047 
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TABLE III. Sensitivities achievable at the 95% confidence level for the anomalous ZZ^ cou- 
plings /if and hf in a) pp — ► Zj + X — > ^ + £~7 + X and b) pp — > Z7 + X — > ^ T 7 + X, at the LHC 
(i/5 = 14 TeV). Results are shown for integrated luminosities of 10 fb _1 and 100 fb _1 . The limits 
for each coupling apply for arbitrary values of the other coupling listed in this table. The CP 
violating couplings hf 2 and all Z77 couplings are assumed to vanish. For the form factor we use 
the form of Eq. (|j) with n = 3 (n = 4) for hf (h z ) and App = 3 TeV. The cuts summarized in 
Sec. IIIB are imposed. For the jet definition, we have used Eq. (|l~0|). 





a) - 


-> Z7 + X -» £+^~ 7 + X 




coupling 


LO 


JXdi = 10 fb- 1 

NLO incl. 


NLO 0-jet 


\h z 1 
l n 30l 


4.6 x 10~ 3 


5.0 x 10~ 3 


4.7 x 10~ 3 


\h z 1 


3.6 x 10~ 5 


3.9 x 10~ 5 


3.7 x 10~ 5 


coupling 


LO 


JCdt = 100 fb _1 

NLO incl. 


NLO 0-jet 


\h z 1 


2.5 x 10~ 3 


2.8 x 10~ 3 


2.6 x 10~ 3 


\h z 1 


1.7 x 10~ 5 


1.9 x 10~ 5 


1.8 x 10~ 5 




b) 


Z 1 + x ^ i> Tl + X 




coupling 


LO 


JCdt = 10 fb" 1 

NLO incl. 


NLO 0-jet 


\h z 1 


3.4 x 10~ 3 


3.7 x 10- 3 


3.5 x 10~ 3 


\h z 1 


2.5 x 10~ 5 


2.7 x 10~ 5 


2.6 x 10~ 5 


coupling 


LO 


JCdt = 100 fb" 1 

NLO incl. 


NLO 0-jet 


l n 3ol 


1.9 x 10 -3 


2.0 x 10 -3 


1.9 x 10~ 3 


\h z 1 

l rt 40l 


1.2 x 10~ 5 


1.3 x 10~ 5 


1.2 x 10~ 5 
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FIGURES 




FIG. 1. Feynman diagrams for the Born level process qq — ► Z7 in the Standard Model. 
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Y 



FIG. 2. Additional Feynman diagrams which contribute to the Born level process qq — ► Z7 
when non-standard model ZZ7 and Z77 couplings are introduced. 
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FIG. 3. Feynman rule for the general ZjV (V = Z, 7) vertex. The factor e is the charge of the 
proton. The vertex function T^y(qi, (fe, P) is given in Eq. (|l|). 
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FIG. 4. Differential cross sections versus Pxil) f° r ( a ) PP —> %1 + X — > i + t 7 + X at 
y/s = 1.8 TeV, and (b) pp -> Z7 + X -> l+i'j + X at yfe = 14 TeV in the SM. The jet-inclusive 
cross sections are shown at the Born level (dashed curves) and with the NLO corrections (solid 
curves). The cuts imposed are summarized in Sec. IIIB. 
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FIG. 5. The differential cross section for the photon transverse momentum in the reaction 
pp — > Z7 + X — > £ + £~j + X at -y/s = 1.8 TeV, (a) in the Born approximation and (b) including 
NLO QCD corrections. The curves are for the SM (solid), /if = 1.0 (dashed), and hf = 0.05 
(dotted). The cuts imposed are summarized in Sec. IIIB. 
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FIG. 6. The differential cross section for the photon transverse momentum in the reaction 
pp — ► Z7 + X — > + X at y/s = 14 TeV, (a) in the Born approximation and (b) including 

NLO QCD corrections. The curves are for the SM (solid), /if = 0.01 (dashed), and hf = 1 x 10~ 4 
(dotted). The cuts imposed are summarized in Sec. IIIB. 
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FIG. 7. Ratio of the NLO and LO differential cross sections of the transverse momentum of 
the photon as a function of pril) for (a) pp — > + X — ► £ + £ _r y + X at ^/s = 1.8 TeV, and (b) 

-> Z7 + X -> £+^7 + X at v/s = 14 TeV. The solid curves show the SM result. The dashed 
and dotted lines display the cross section ratio for non-zero values of /if and hf , respectively. The 
cuts imposed are summarized in Sec. IIIB. 
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FIG. 8. The pt{i) differential cross section for (a) pp — > Zj+X — > £ + £ ^y+X at ^/s = 1.8 TeV, 
and (b) pp Z7 + X -> ^+^"7 + X at ^/s = 14 TeV in the SM. The inclusive NLO differential 
cross section (solid line) is decomposed into the 0(a s ) 0-jet (dotted line) and LO 1-jet (dashed 
line) exclusive differential cross sections. For comparison, the Born cross section (dot dashed line) 
is also shown. The cuts imposed are summarized in Sec. IIIB. For the jet definitions, we have used 



Eqs. © and (g(J). 
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FIG. 9. The pr(l) differential cross section for (a) pp — * Zj + jet — » t + l~~{ + jet at 
yfs = 1.8 TeV, and (b) pp — > Z7 + jet — > ^ + £~7 + jet at -y/s = 14 TeV. The curves in part (a) 
are for the SM (solid), h§ = 1.0 (dashed), and hf = 0.05 (dotted). In part (b), the dashed and 
dotted curves are for h% = 10~ 2 and hf = 10~ 4 , respectively. The cuts imposed are summarized 
in Sec. IIIB. For the jet definitions, we have used Eqs. (^) and (IH). 
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FIG. 10. (a) The pr{l) differential cross section, and (b) the ratio of the NLO 0-jet to 
LO differential cross sections of the photon transverse momentum as a function of Pt(j) for 
pp -> Z7 + jet -> ^ T 7 + jet at y/s = 1.8 TeV. The curves are for the SM (solid), /if = 1.0 
(dashed), and hf = 0.05 (dotted). The cuts imposed are summarized in Sec. IIIB. 
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FIG. 11. (a) The pr{l) differential cross section, and (b) the ratio of the NLO 0-jet to 
LO differential cross sections of the photon transverse momentum as a function of pt(j) for 
pp -> Z7 + jet -> /5 T 7 + jet at sfs = 14 TeV. The curves are for the SM (solid), /if = 1CT 2 
(dashed), and hf = 10 -4 (dotted). The cuts imposed are summarized in Sec. IIIB. 
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FIG. 12. Limit contours at the 95% CL for pp — > Zj-\-X — » £ + £~j + X (I = e, fx), derived from 
the pr(l) distribution at the Tevatron for JCdt = 1 fb _1 . Contours are shown in three planes: a) 
the fofo - hj plane, b) the h§ - hj plane, and c) the hf - hj plane. The solid lines give the 
results for LO Zj production. The dashed curves give the inclusive NLO results, and the dotted 
lines show the bounds obtained from the exclusive Zj + jet channel. In each graph, only those 
couplings which are plotted against each other are assumed to be different from their zero SM 
values. The cuts imposed are summarized in Sec. IIIB. For the jet definition, we have used Eq. @. 
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FIG. 13. Limit contours at the 95% CL in the /if - hf plane for pp -> Z7 + X -> £+£"7 + X 
(I = e, ^u), derived from the pt{i) distribution at the Tevatron for a) jCdt = 1 fb _1 and b) 
/ Cdt = 10 fb . The solid lines give the results for LO Z7 production. The dashed curves give the 
inclusive NLO results, and the dotted lines show the bounds obtained from the exclusive Z7 + O jet 
channel, hf 2 arid all Z77 couplings are assumed to be zero. The cuts imposed are summarized in 
Sec. IIIB. For the jet definition, we have used Eq. (|9|). 
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FIG. 14. Limit contours at the 95% CL in the h% Q — hf plane for pp — > Z7 + X — > j4 T 7 + X 
derived from the pt{i) distribution at the Tevatron for /£<ii = 1 fb _1 . The solid line gives 
the result for LO Z7 production. The dotted curve shows the bounds obtained from the NLO 
calculation. hf 2 and all Z77 couplings are assumed to be zero. The cuts imposed are summarized 
in Sec. IIIB. 
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FIG. 15. Limit contours at the 95% CL in the h§ - hf plane for pp -> Z7 + X -> £+£"7 + X 
(•£ = e, /i) derived from the pr{l) distribution at the LHC for fCdt = 100 fb _1 . The solid line 
gives the results for LO Z7 production. The dashed curve gives the inclusive NLO results, and 
the dotted lines shows the bounds obtained from the exclusive Z7 + jet channel, hf 2 and all 
Z77 couplings are assumed to be zero. The cuts imposed are summarized in Sec. IIIB. For the jet 



definition, we have used Eq. (10). 
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